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Hydrodynamiczne aspekty przepływu cytoplazmy w międzywęzłowych 
komórkach Characeae 


TuaqpoanHamnueckue ACNeKTbI TEYEHHA LUKTONIIA3MbI B MEKĄDYZNEHHbIX KNETKAX Characeae 


Some Hydrodynamical Aspects of Cytoplasmic Streaming in the Imternodal Cells 
of Characeae 


Ruch cytoplazmy w międzywęzłowych komórkach Characeae [1] wy- 
kazuje prawidłowości, dzięki którym możliwe jest przeprowadzenie pew- 
nych obliczeń opartych na hydrodynamice cieczy lepkich. W komórkach 
tych przepływ cytoplazmy (Ściślej, wewnętrznej jej warstwy zwanej en- 
doplazmą) i soku komórkowego zachodzi wzdłuż dwu zwiniętych ze so- 
bą półwalców, oddzielonych pasmem ciszy, w którym prędkość przepły- 
wu jest równa zeru (ryc. 1). W obu półwalcach prędkości przepływu ma- 
ją kierunki przeciwne. 

Endoplazma i sok komórkowy są cieczami o pewnej lepkości, aby więc 
raogły znajdować się w ciągłym ruchu, muszą na nie działać stale jakieś 
siły. Istnieje wiele teorii wyjaśniających w sposób jakościowy przepływ 
endoplazmy (1, 4, 6); obecnie przyjęty jest pogląd, według którego ,,sila 
napędowa” umiejscowiona jest w pobliżu ścianki komórkowej, na grani- 
cy dwu warstw cytoplazmy: nieruchomej (ektoplazmy) i płynącej (endo-- 
plazmy), stycznie do ich powierzchni. Za takim umiejscowieniem siły na- 
pędowej przemawia rozkład prędkości przepływu w różnych punktach 
komórki; prędkość ta osiąga największą wartość w pobliżu ścianki komór- 
kowej. Tazawa i Kishimoto [2] zmierzyli wartość siły napedo- 
wej, doprowadzając endoplazme do zatrzymania się, wskutek przepływu 
cieczy wzdłuż wodniczki uciętego na końcach międzywęźla. Wyznaczona 
w ten sposób wartość siły napędowej (równej w tym wypadku sile stycz- 
nej do endoplazmy i zatrzymującej jej ruch) wynosiła 1—2 dyn/cm?. K a- 
miya i Kuroda [3] otrzymali taką samą wartość siły napędowej, do- 
prowadzając endoplazmę do zatrzymania w mikroskopie wirówkowym. 
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W niniejszej pracy przeprowadzono obliczenia mające na celu znale- 
zienie ilościowej zależności między siłą napędową, lepkością endoplazmy 
i prędkością jej ruchu w różnych punktach komórki. Do przeprowadze- 
nia obliczeń konieczne było wprowadzenie uproszczenia, polegającego na 
zaniedbaniu śrubowo-spiralnego charakteru przepływu i na przyjęciu, że 
zachodzi on wzdłuż dwu półwalców nie zwiniętych ze sobą, lecz tylko 
równolegle do siebie przylegających (ryc. 2). Ze względu na to, że skok 
spirali, jaką tworzą zwinięte ze sobą półwalce (patrz ryc. 1) jest bardzo 


pas ciszy 


cytoplazma 


sok komórkowy f 


Ryc. 1. Schemat budowy międzywęźla glonu Characeae z zaznaczonym przekrojem 
poprzecznym; strzałki wskazują kierunek przepływu endoplazmy i soku komórkowego 


duży w porównaniu ze średnicą komórki, powyższe założenie nie ma prak- 
tycznie żadnego wpływu na wyniki ilościowe obliczeń, upraszcza nato- 
miast znacznie matematyczną stronę zagadnienia. 


ROWNANIE RUCHU ENDOPLAZMY I SOKU KOMÓRKOWEGO 


Niech w cylindrze o promieniu R (patrz ryc. 2) znajdują się dwie cie- 
cze nie mieszające się ze sobą (sok komórkowy i endoplazma) o współ- 
czynnikach lepkości odpowiednio nı i n». Jedna z tych cieczy wypełnia 
wewnętrzny walec o promieniu r, (sok komórkowy), druga obszar pomię- 
dzy tym walcem a walcem zewnętrznym o promieniu R (endoplazma). Za- 


Ryc. 2. Uproszczony schemat międzywęźla z zaniedbaniem spiralnego charakteru 

ruchu endoplazmy; zaznaczono współrzędne walcowe z, r, œ i współczynniki lepkości 

Y i qe; R i Tę — odpowiednio zewnętrzny i wewnętrzny promień warstwy endo- 
plazmy 
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kładamy, że obie ciecze są nieściśliwe i ciśnienie wewnątrz nich jest wszę- 
dzie jednakowe. Siły poruszające obie ciecze niech działają stycznie do 
zewnętrznej powierzchni endoplazmy (strzałki na ryc. 2). 

Przy rozwiązywaniu równania ruchu tych cieczy użyjemy współrzęd- 
nych walcowych (oś Z jest skierowana wzdłuż osi walca, kąt o jest kątem 
pomiędzy pasem ciszy, gdzie prędkość przepływu jest równa zeru, a kie- 
runkiem promienia wodzącego r). Ogólne równanie ruchu cieczy lepkiej, 
nieściśliwej (równanie Naviera — Stokesa) ma postać [7]: 


pat +AVV)V -nAV = -grad p 
(gdzie: go — gęstość cieczy, V — wektor prędkości, 17 — współczynnik lep- 
kości, p — ciśnienie). 
Przy stacjonarnym przepływie cieczy jej prędkość nie zależy od cza- 


su ani od współrzędnej z. W tych warunkach składowe radialna i kąto- 
wa wektora prędkości są równe zeru: 


V=[V-,Ve,Ve] =[o,0virg)] 


Tak więc równanie Naviera — Stokesa Pera kształt następujący: 


PRV + rze + =0 (1) 


Przeprowadzamy rozdzielenie wc przedstawiając poszukiwa- 
ną funkcję V(r, Q) w postaci iloczynu: 
V(r.) Plr) PQ) 
Otrzymujemy w ten sposób dwa równania: 


p Er AO 
dr? T" dr F 


aod 
dą? = a? © 
gdzie a jest stałą rozdzielenia zmiennych. 
Ogólnym rozwiązaniem równania (1) jest funkcja: 
V (r, g)=(Ar*+ Br"*) (C sin a ọ+ D cos a q). 
Funkcja V (r, Q) musi spełniać następujące warunki symetrii (patrz ryc. 2): 
V (r, ptx)=—V (r, 9) 
V(r, —9)=—V A G?) 
V (r, 0)=V (r, x)=0 
Wynika stąd, że: 


4 Annales. sectio AA, vol. XXVIII 
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u/143, 5505 leb; a. 


D=0 
Na funkcję V (r, Q) mamy więc nieskończenie wiele rozwiązań typu: 
V (r, 9)=(A1 r” +B’, r”) sin nę (2) 


gdzie n oznacza kolejne liczby nieparzyste. Ponadto w obszarze wewnętrz- 
nego walca zajętego przez ciecz o lepkości nı (sok komórkowy, patrz ryc. 
2) występuje punkt o współrzędnej r=0, a zatem w rozwiązaniu odpowia- 
dającym temu obszarowi musi być B,=0. 

Tak więc dla cieczy wypełniającej walec zewnętrzny o przekroju pier- 
ścieniowym (endoplazma) mamy rozwiązanie ogólne w postaci sumy roz- 
wiązań danych wzorem (2): 


V4\r.9) =), (Anr?+ Bqr""lsin np 


n=1.3,5.. 
zaś dla cieczy wypełniającej wewnętrzny walec: 


virg) EY Car?sin np 
n «43.5... 
Dla udogodnienia dalszych rachunków podstawiamy: 
pal 
R De | 


i otrzymujemy: 


=x , =x% , ApRT=An: BhR=Bq: CnR=Cn 


vlx.ę) = Teru ng 


n=4,3,5... g 
(3) 
K(x9) = S lan + Bn xo) sin np 

nai 3.5... 
Korzystając z ogólnego rozwiązania równania ruchu dwu cieczy w prze- 
wodzie walcowym danego wzorami (3), rozpatrzymy trzy warianty tego 
ruchu różniące się odmiennymi założeniami dotyczącymi mechanizmu 
napędowego. 


Wariant I (stała prędkość) 


Zakładamy, że na powierzchni płynącej endoplazmy działa taki me- 
chanizm napędowy, który zmusza cząstki endoplazmy do ruchu ze sta- 
łą prędkością, niezależną od kąta Q, tj. od odległości od pasa ciszy. 
Wszystkie cząsteczki na zewnętrznej powierzchni górnego półwalca (ryc. 
2) posiadają prędkość +V zaś na powierzchni dolnego półwalca — Vo. 
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Przyjęte w tym wariancie założenie o stałej prędkości na powierzchni 
endoplazmy narzuca następujący warunek brzegowy na wzory (3): 


yda OX P<KT 
x=1,9) = (4) 
a - dla <p [L20 
W odległości ry od osi komórki (powierzchnia wewnętrznego walca, ryc. 
2), czyli dla x=2,) sok komórkowy i endoplazma oddziaływają na siebie 


jednakowymi siłami i prędkości ich są jednakowe. Warunek równości sił 
zapisujemy w postaci: 
=p 
Ró 
c 


ER 


nee (5) 


ER 


przy obliczaniu pochodnych prędkości względem r posługujemy się za- 
leżnością: 

av _av ox _1.9V 

Arin el! Opec ike va 


Warunek równości prędkości: 
V (Xo, 9)= Və (Xo, q) (6) 


Podstawiając do zależności (4), (5) i (6) wzory (3) otrzymujemy następują- 
ce równania: 


veering da 0Xo<n (4a) 
n=4.2.5... =K dla TU <9 <zn 
MEn" = Tal An x0 = By x0") 

(5a) 


n n -n 
CnX = An%+ Bn x; 


Z równania (4a) widać, że wyrażenie (A”+ B„) musi być równe n-temu 
współczynnikowi rozwinięcia w szereg Fouriera funkcji: 


Vig) = %dla 0<o< T 
Wda r<p<2n 


czyli, że 
4V 
A = 0 4b 
n+ By nn (4b) 
Rozwiązując układ równań (4b), (5a) i (6a) ze względu na współczyniki 
An, Bn, Cn, otrzymujemy po ich podstawieniu do zależności (3) następują- 
ce wzory na prędkości soku komórkowego i endoplazmy: 
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JĄ (1+5)x7 ©. 
V(x,9) 5% dż 2 niren sin NQ 
h ei. x". 
Klx,9) = 2K y walc sin np (7) 
n=4,2,5,... 
gdzie: ô= Ee 


Wariant II (stała siła napędowa) 


Załóżmy, że na zewnętrzną powierzchnię endoplazmy działa stała 
siła napędowa Fo, przy czym jej zwrot jest wszędzie zgodny ze 
zwrotem prędkości. Mówiąc o stałości siły napędowej, mamy na myśli jej 
niezależność od kata Q w przedziale jego wartości od zera do x oraz od x 
do 2x (siła, podobnie jak w poprzednim wariancie prędkość, w obszarze 
pasa ciszy zmienia kierunek na przeciwny). 

Warunek brzegowy dla prędkości V2 (x, Q) na powierzchni endoplazmy 
ma wtedy postać następującą: 


-Adlan <9 <27 


Pozostałe dwa warunki brzegowe są takie same jak w poprzednim wa- 
riancie (równania (5a) i (6a)). Równanie (8) po podstawieniu doń zależ- 
ności (3) przybiera postać: 


(8) 


= F 
) ml — 284] sin ng = g dia 0<g<n (8a) 
n=4,3,5,. R ai 74 dla nrKp<2nr 


Stąd widać, że wyrażenie R E - ca musi być równe n-temu współ- 


czynnikowi rozwinięcia w szereg Fouriera funkcji: 


KR da Ko <N 


SAAB (e dla rKo<2n 


Tn 


a więc: h [240 = a MAC (8b) 


Rozwiązując układ równań (8b), (5a) i (6a), otrzymujemy następujące 
wzory na V; (x, g)i V2 (x, Q): 


Hydrodynamiczne aspekty przepływu cytoplazmy... 53 
LRR (1+6) x 
K(x.p) = Nh PE (1-82 sin np 
x n 
Valx,p) = iga Sasa sin ng 


gdzie: 6 


(9) 


najn A m". 


Wariant III (siła napędowa zależna od prędkości) 


Przyjmujemy, że siła napędowa działająca w dowolnie wybranym 
punkcie powierzchni endoplazmy zależy od prędkości, jaką en- 
doplazma posiada w tym punkcie. 

Maksymalną prędkość osiąga endoplazma w punktach o współrzęd- 
nych x=1, p= +x/2. Oznaczmy ją przez Vo. Załóżmy, że. siła napędowa 
działająca na endoplazmę w danym punkcie jej powierzchni jest tym 
większa, im bardziej różni się prędkość endoplazmy w tym punkcie od 
V,. Ilościowe ujęcie tej zależności przedstawia następujące równanie: 


F(R,g) =k[V-VIRgl] (10) 


Jak widać z tego równania, siła napędowa powinna być najmniejsza wów- 
czas gdy V (R, g)=V, (to jest w punktach o współrzędnych x=1, p= +21/2) 
zaś największa w obrębie pasa ciszy i w jego bezpośrednim sąsiedztwie, 
wtedy bowiem V (R, ọ) jest najbardziej zbliżone do zera. Podobnie w cza- 
sie regeneracji ruchu strumienia endoplazmy, np. po pobudzeniu komór- 
ki, siła napędowa powinna być największa w początkowej fazie procesu 
regeneracji, gdy prędkość endoplazmy jest jeszcze mała, zaś zbliża się 
do zera, gdy prędkość endoplazmy osiąga swą maksymalną wartość Vo. 
Taki charakter siły napędowej wydaje się wyjaśniać, dlaczego przyspie- 
szenie ruchu endoplazmy po pobudzeniu komórki jest największe w po- 
czątkowym etapie procesu regeneracji, zaś maleje w miarę, jak prędkość 
endoplazmy osiąga wartość V, [1]. 

Dopasujmy ogólne rozwiązanie równania ruchu (3) do nowych warun- 
ków brzegowych. Na powierzchni endoplazmy (x=1) działa siła styczna 
równa: 


(11) 


parka dla 0< <v 


F 3 = 2 
is 4 -ĘG-kVlRg) dia n < <21 


r=R 
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czyli: 
7. Br). R-k ElAp+Bh) sin ny 
772 [ 242 - 7 |sinnę = 7 
Stąd zaś: 
= K dla O <g<n 
nAn _ aBa) + klA +B, |sinn = o 
l PRA WPW |e dla n<ę<2n 


Podobnie jak w poprzednim wariancie, uzyskujemy równanie: 


n( 240 og) Hanan = 48 


Rozwiązując układ równań (12), (5a) i (6a), otrzymujemy na V, (a, @) oraz 
V2 (x, p) następujące wzory: - 


(12) 


—_ LRRII+6 E x” : 
Valxig) = Abp ) FAIR RAGE ny PK 
N= 1,3,5,... i 
| = 4GR A z RZE aa a has 
KU T r l1-0%")+kRnl1i+5% |) 


N Ya 
$- 


(13) 


POROWNANIE DANYCH TEORETYCZNYCH Z DOSWIADCZALNYMI 


Wyniki przedstawionych powyżej obliczeń porównano z doświadczal- 
nymi krzywymi rozkładu prędkości w komórce międzywęzłowej Chara 
brauni. Używane do doświadczeń komórki miały średnicę 0,4 mm i dłu- 
gość 3,5—5 cm. Prędkość przepływu endoplazmy i soku komórkowego 
mierzono posługując się mikroskopem i stoperem. Pomiary przeprowadzo- 
no śledząc ruch możliwie najmniejszych cząstek zawieszonych w płynącej 
endoplazmie. 

Określono radialny rozkład prędkości endoplazmy i soku komórkowe- 
go dla Q=au/2, to jest w płaszczyźnie leżącej w połowie odległości między 
pasami ciszy (ryc. 3a) oraz kątowy rozkład prędkości w endoplazmie 
(ryc. 3b). Przy wyznaczaniu kątowego rozkładu prędkości ustawiano stre- 
ię głębi ostrości mikroskopu tuż pod ścianką komórki tak, aby mierzyć 
prędkości przepływu jak najbardziej zewnętrznych warstw endoplazmy. 
Kąt @ wyliczano ze wzoru sin Q=d/R, gdzie R — promień komórki, zaś d — 
odległość obserwowanego punktu od pasa ciszy (patrzy ryc. 3b). Oblicza- 
jąc teoretyczne wartości prędkości przepływu, przyjęto, że grubość war- 
stwy endoplazmy wynosi 15 u. 
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Ryc. 3. a — schemat pomiaru radialnego rozkładu prędkości; b — pomiar kątowego 
rozkładu prędkości; M — obiektyw mikroskopu, CC’ — pas ciszy, B, ... Bn — punkty, 
w których mierzono prędkości przepływu 


Dla promienia powierzchni R=200 u otrzymujemy wtedy wartość 
promienia wodniczki r,= 185 u, co daje x= 0,925. 

Do wzorów na prędkości przepływu otrzymanych w wariancie I współ- 
czynnik lepkości endoplazmy wchodzi tylko za pośrednictwem parame- 
tru 0, zaś we wzorach wariantów II i III występuje on także w sposób 
jawny. Przy obliczeniach teoretycznych wartość współczynnika lepkości 
endoplazmy dobierano tak, aby w każdym wariancie przy wartości siły 
napędowej w granicach 1—2 dyn/cm? [1, 2, 3] teoretyczne wartości V,=V 
(1, 1/2) zgadzały się z doświadczalnie znalezioną wartością V,=70 u/sek. 
Założono ponadto, że ’=1, co przy stosowanych w tej pracy wartościach 
nı i n2 (M1 = 0,01 P, n: > 0,5 P) nie wprowadza do wartości prędkości 
błędu większego niż 1—2%. 

Występujący we wzorach wariantu III parametr k jest możliwy do 
wyznaczenia na podstawie danych uzyskanych w eksperymentach prze- 
prowadzonych przez Tazawę i Kishimoto [2] oraz Kamiya 
i Kuroda [3]. W eksperymentach tych mierzona była siła napędowa 
występująca wtedy, gdy endoplazma jest zatrzymana siłami zewnętrzny- 
mi. Mierzona w takich warunkach siła, w świetle przyjętych tu założeń 
(10) jest maksymalną wartością siły, z jaką mechanizm napędowy zdol- 
ny jest oddziaływać na endoplazmę. W przypadku V (R,q)=0 (co odpo- 
wiada warunkom tych doświadczeń) równanie (10) przyjmuje postać: 


F (R, q)=kV,=F, 
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10 20 30 40 so 60 70 V (u/sek.) 


Ryc. 4. Radialne rozkłady prędkości; 1, 2, 3 — krzywe teoretyczne, 4 — punkty doś- 
wiadczalne; zaznaczono wartości ne przyjęte przy obliczaniu krzywych wariantu III 


V (u/sek! 
80 
—— M,05P 


70 0— [i,1P 


60 
50 
40. 
30 
20 


10 


e ee OT RTZICZ 
o KA 10° 20° 309 40° 50° 60° 70° 80° 90° Y 
ciszy (2°) 
Ryc. 5. Kątowe rozkłady prędkości; 1 — punkty doświadczalne; zaznaczono wartości 
Ng użyte przy obliczaniu krzywych wariantu III 
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Stąd otrzymujemy: 


Tazawa i Kishimoto przyjmowali wartość Fo/Vo za miarę lepkości; wyz- 
naczona z uzyskanych przez nich danych średnia wartość k wynosi 
214 g/cm? sek., przy średniej wartości F,= 1,5 dyn/cm?. 

Porównanie wyników doświadczalnych z teoretycznymi przedstawia- 
ją ryc. 4 i 5. Jak widać z ryc. 4, radialne rozkłady prędkosci uzyskane 
w wariancie I i III mechanizmu napędowego dosyć dobrze zgadzają się 
z punktami doświadczalnymi; krzywa uzyskana w wariancie II nieco od 
nich odbiega. Na ryc. 5 widoczna jest stosunkowo dobra zgodność krzy- 
wej wariantu III z danymi doświadczalnymi, natomiast krzywe teoretycz- 
ne wariantu I i II całkowicie odbiegają od punktów eksperymentalnych. 
Postulowane w wariancie I własności mechanizmu napędowego w obsza- 
rze małych kątów (tj. w pobliżu pasa ciszy) są nieco wyidealizowane, to- 
też pewnej niezgodności wariantu I z doświadczeniem w tym obszarze 
należało się spodziewać. Jednak, jak widać z ryc. 5, rozbieżności między 
krzywą odpowiadającą wariantówi I a punktami doświadczalnymi rozcią- 
gają się w obszarze około siedmiokrotnie szerszym od pasa ciszy. 

O nierealności wariantu II (założenie stałej siły napędowej), oprócz 
rozbieżności widocznych na ryc. 4 i 5, świadczy też fakt, że w wariancie 
tym dla uzyskania prawidłowej wartości maksymalnej prędkości: V na- 
leżało użyć współczynnika lepkości endoplazmy równego 40—80 P (przy 
wartości siły napędowej 1—2 dyn/cm?). Wartości n» cytowane w piśmien- 
nictwie wahają się w granicach od ułamka puaza do kilku puazów [5]. 


PODSUMOWANIE 


W świetle otrzymanych wyników najsłuszniejszy wydaje się taki me- 
chanizm napędowy ruchu cytoplazmy, który wywiera na endoplazmę siłę 
zależną od prędkości, jaką endoplazma posiada (równanie (10)). 

Zależność siły napędowej od prędkości przepływu endoplazmy stawia 
w nowym świetle doświadczenia Tazawy i Kishimoto [2] oraz 
Kamiya i Kuroda [3], mające na celu wyznaczenie wartości siły 
napędowej. Siła ta może być równa 1—2 dyn/cm?, jeśli endoplazma zo- 
staje doprowadzona do zatrzymania się siłami zewnętrznymi. Natomiast 
— zgodnie z założeniami wariantu III — może mieć ona wartość zupełnie 
inną (mniejszą), gdy endoplazma znajduje się w normalnym ruchu. 

Związek parametru k z innymi wielkościami charakteryzującymi prze- 
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pływ cytoplazmy pozostaje nie wyjaśniony. Wynika to z fenomenologicz- 
nego charakteru przeprowadzonych rozważań, nie wnikających w mikro- 
skopową naturę mechanizmu napędowego. 
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PE3IOME 


B pa6oTe nonyueHbi ypaBH€HMA CBAJbIBAIOLIME ABWKYLYIO CHNY NPOTO- 
nna3Mbl, CKOPOCTŁ TEUEHMA H BA3KOCTb. CPaBHEHKE TEOpPETMUECKHX n 3KCNE- 
PMUMEHTANbHbIX MAHHbIX YKASbIBACT Ha 3ABHCHMOCTŁ ĄBUKYLLEŃ CHNbI OT 
CKOPOCTM TEHUEHKA' 


SUMMARY 


In the presented paper some equations were obtained, which con- 
tain the motive force of protoplasmic streaming, velocity of streaming 
and viscosity of protoplasm. The comparison of the theoretical and expe- 
rimental data indicates that the motive force of protoplasmic streaming 
depends on the velocity of streaming. 


